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A STOCHASTIC-LAGRANGIAN APPROACH TO THE 
NAVIER STOKES EQUATIONS IN DOMAINS WITH BOUNDARY 

By Peter Constantin^ and Gautam Iyer^ 

University of Chicago and Carnegie Mellon University 

In this paper we derive a probabilistic representation of the de- 
terministic 3-dimensional Navier-Stokes equations in the presence of 
spatial boundaries. The formulation in the absence of spatial bound- 
aries was done by the authors in [Comm. Pure Appl. Math. 61 (2008) 
330-345]. While the formulation in the presence of boundaries is sim- 
ilar in spirit, the proof is somewhat different. One aspect highlighted 
by the formulation in the presence of boundaries is the nonlocal, im- 
plicit influence of the boundary vorticity on the interior fluid velocity. 

1. Introduction. The (unforced) incompressible Navier-Stokes equations 

(1.1) dtu+ {u-V)u- uAu + Vp = 0, 

(1.2) V-u = 

describe the evolution of the velocity field u of an incompressible fluid with 
kinematic viscosity > in the absence of external forcing. Here u = u{x, t) 
with t>0, X ^W^, d>2. Equation (1.2) is the incompressibility constraint. 
Unlike compressible fluids, the pressure p in (1.1) does not have a physical 
meaning and is only a Lagrange multiplier that ensures incompressibility is 
preserved. While equations (1.1) and (1.2) can be formulated in any dimen- 
sion d>2, they are usually only studied in the physically relevant dimen- 
sions 2 or 3. The presentation of the Navier-Stokes equations above is in 
the absence of spatial boundaries; an issue that will be discussed in detail 
later. 
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When 1/ = 0, (1.1) and (1.2) are known as the Euler equations. These 
describe the evolution of the velocity field of an (ideal) inviscid and incom- 
pressible fluid. Formally the difference between the Euler and Navier-Stokes 
equations is only the dissipative Laplacian term. Since the Laplacian is ex- 
actly the generator a Brownian motion, one would expect to have an exact 
stochastic representation of (1.1) and (1.2) which is physically meaningful, 
that is, can be thought of as an appropriate average of the inviscid dynamics 
and Brownian motion. 

The difficulty, however, in obtaining such a representation is because of 
both the nonlinearity and the nonlocality of equations (1.1) and (1.2). In 2D, 
an exact stochastic representation of (1.1) and (1.2) dates back to Chorin [14] 
in 1973 and was obtained using vorticity transport and the Kolmogorov 
equations. In three dimensions, however, this method fails to provide an 
exact representation because of the vortex stretching term. 

In 3D, a variety of techniques has been used to provide exact stochastic 
representations of (1.1) and (1.2). One such technique (Le Jan and Sznit- 
man [26]) uses a backward branching process in Fourier space. This ap- 
proach has been extensively studied and generalized [3, 4, 32, 35, 36] by 
many authors (see also [37]). A different and more recent technique due 
to Busnello, Flandoli and Romito [6] (see also [5]) uses noisy flow paths 
and a Girsanov transformation. A related approach in [11] is the stochastic- 
Lagrangian formulation, exact stochastic representation of solutions to (1.1) 
and (1.2) which is essentially the averaging of noisy particle trajectories 
and the inviscid dynamics. Stochastic variational approaches (generalizing 
Arnold's [1] deterministic variational formulation for the Euler equations) 
have been used by [13, 16] and a related approach using stochastic differen- 
tial geometry can be found in [19]. 

One common setback in all the above methods is the inability to deal 
with boundary conditions. The main contribution of this paper adapts the 
stochastic-Lagrangian formulation in [11] (where the authors only consid- 
ered periodic boundary conditions or decay at infinity) to the situation with 
boundaries. The usual probabilistic techniques used to transition to domains 
with boundary involve stopping the processes at the boundary. This intro- 
duces two major problems with the techniques in [11]. First, stopping intro- 
duces spatial discontinuities making the proof used in [11] fail and a different 
approach is required. Second and more interesting is the fact that merely 
stopping does not give the no-slip (0-Dirichlet) boundary condition as one 
would expect. One needs to also create trajectories at the boundary which 
essentially propagate the influence of the vorticity at the boundary to the 
interior fluid velocity. 

1.1. Plan of the paper. This paper is organized as follows. In Section 2 
a brief introduction to the stochastic-Lagrangian formulation without bound- 
aries is given. In Section 3 we motivate and state the stochastic-Lagrangian 
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formulation in the presence of boundaries (Theorem 3.1). In Section 4 we 
recall certain standard facts about backward Ito integrals which will be used 
in the proof of Theorem 3.1. In Section 5 we prove Theorem 3.1. Finally, in 
Section 6 we discuss stochastic analogues of vorticity transport and inviscid 
conservation laws. 

2. The stochastic-Lagrangian formulation without boundaries. In this 
section, we provide a brief description of the stochastic-Lagrangian formu- 
lation in the absence of boundaries. For motivation, let us first study a La- 
grangian description of the Euler equations [equations (1.1) and (1.2) with 

= 0; we will usually use a superscript of to denote quantities relating to 
the Euler equations]. Let d = 2,3 denote the spatial dimension and be 
the flow defined by 

(2.1) x^=unxh 

with initial data ^o('^) ~ ^ ^ clarify our notation, X'^ is 

a function of the initial data a € M*^ and time t € [0, oo). We usually omit the 
spatial variable and use X^ to denote X^{-,t), the slice of X^ at time t. Time 
derivatives will always be denoted by a dot or dt instead of a t subscript. 

One can immediately check (see, e.g., [7]) that u satisfies the incompress- 
ible Euler equations if and only if X^ is a gradient composed with X. By 
Newton's second law, this admits the physical interpretation that the Euler 
equations are equivalent to assuming that the force on individual particles 
is a gradient. 

One would naturally expect that solutions to the Navier-Stokes equations 
can be obtained similarly by adding noise to particle trajectories and aver- 
aging. However, for noisy trajectories, an assumption on X^ will be prob- 
lematic. In the incompressible case, we can circumvent this difficultly using 
the Weber formula [38] [equation (2.2) below]. Indeed, a direct computation 
(see, e.g., [7]) shows that for divergence free u, the assumption that X^ is 
a gradient is equivalent to 

(2.2) nO = P[(VMO)(nOoAO)], 

where P denotes the Leray-Hodge projection [10, 15, 28] onto divergence 
free vector fields, the notation V* denotes the transpose of the Jacobian 
and for any t > 0, = {X^ )~^ is the spatial inverse of the map X^ [i.e., 
A°{Xf{a)) = a for ah a G and X^{Af{x)) = x for all x € M"^]. 

From this we see that the Euler equations are formally equivalent to equa- 
tions (2.1) and (2.2). Since these equations no longer involve second (time) 
derivatives of the fiow X^, one can consider noisy particle trajectories with- 
out any analytical difficulties. In fact, adding noise to (2.1) and averaging 
out the noise in (2.2) gives the equivalent formulation of the Navier-Stokes 
equations stated below. 
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Theorem 2.1 (Constantin, Iyer [11]). Let d£ {2,3} be the spatial di- 
mension, u > represent the kinematic viscosity and uq be a divergence free, 
periodic, Holder 2 + a function and W be a d-dimensional Wiener process. 
Consider the system 

(2.3) dXt = ut{Xt) dt + V2^dWt, 

(2.4) Xo{a)=a Va G M"^, 

(2.5) ut = EP[{V*At){uooAt)], 

where, as before, for any t>0, At = X^^ denotes the spatial inverse" of Xt- 
Then u is a classical solution of the Navier-Stokes equations (1.1) and (1.2) 
with initial data uq and periodic boundary conditions if and only if u is a fixed 
point of the system (2.3)-(2.5). 

Remark. The flows X,A above are now a function of the initial data 
a G M^, time t G [0,oo) and the probability variable zu G 17. We always sup- 
press the probability variable, use Xt to denote X{-,t) and omit the spatial 
variable when unnecessary. The function ?x is a deterministic function of 
space and time and, as above, we use ut to denote the function u{-,t). 

We now briefly explain the idea behind the proof of Theorem 2.1 given 
in [11] and explain why this method can not be used in the presence of 
spatial boundaries. Consider first the solution of the SDE (2.3) with initial 
data (2.4). Using the Ito-Wentzel formula [25], Theorem 4.4.5, one can show 
that any (spatially regular) process 9 which is constant along trajectories 
of X satisfies the SPDE 

(2.6) dOt + (ut ■ V)et dt - i^AOt dt + V2^Vet dWt = 0. 

Since the process A (which, as before, is defined to be the spatial inverse 
of X) is constant along trajectories of X, the process 9 defined by 

(2.7) 9t = 9oo At 

is constant along trajectories of X. Thus, if is regular enough (C^), 
then 9 satisfies SPDE (2.6). Now, if u is deterministic, taking expected 
values of (2.6) we see that 9t = E9q o At satisfies 

(2.8) dt9t + {ut ■ V)9t - vA9t = 
with initial condition 9\t=Q = 9q. 



^It is well known (see, e.g., Kunita [25]) that the solution to (2.3) and (2.4) gives 
a stochastic flow of difTeomorphisms and, in particular, guarantees the existence of the 
spatial inverse of X. 
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Remark. Note that when u = 0, A is deterministic so 6 = E9 = 9. Fur- 
ther, equation (2.6) reduces to the transport equation for which writing the 
solution as 6t = 9qo At is exactly the method of characteristics. When i/ > 0, 
the above procedure is an elegant generalization, termed as the "method of 
random characteristics" (see [11, 20, 33] for further information). 

Once explicit equations for A and uqo A have been established, a direct 
computation using Ito's formula shows that u given by (2.5) satisfies the 
Navier-Stokes equations (1.1) and (1.2). This was the proof used in [11]. 

Remark. This point of view also yields a natural understanding of gen- 
eralized relative entropies [8, 12, 29, 30]. Eyink's recent work [17] adapted 
this framework to magnetohydrodynamics and related equations by using 
the analogous Weber formula [24, 34]. We also mention that Zhang [39] con- 
sidered a backward analogue and provided short elegant proofs to classical 
existence results to (1.1) and (1.2). 

3. The formulation for domains with boundary. In this section we de- 
scribe how (2.3)-(2.5) can be reformulated in the presence of boundaries. We 
begin by describing the difficulty in using the techniques from [11] described 
in Section 2. 

Let D cR'^ he a domain with Lipschitz boundary. Even if we insist n = 
on the boundary of D, we note that the noise in (2.3) is independent of space 
and thus, insensitive to the presence of the boundary. Consequently, some 
trajectories of the stochastic flow X will leave the domain D and for any 
t > 0, the map Xt will (surely) not be spatially invertible. This renders (2.7) 
meaningless. 

In the absence of spatial boundaries, equation (2.7) dictates that 0{x,t) 
is determined by averaging the initial data over all trajectories of X which 
reach x at time t. In the presence of boundaries, one must additionally av- 
erage the boundary value of all trajectories reaching {x,t), starting on dD 
at any intermediate time (Figure 1). As we will see later, this means the 
analogue of (2.7) in the presence of spatial boundaries is a spatially discon- 
tinuous process. This renders (2.6) meaningless, giving a second obstruction 
to using the methods of [11] in the presence of boundaries. 

While the method of random characteristics has the above inherent dif- 
ficulties in the presence of spatial boundaries, equation (2.8) is exactly the 
Kolmogorov Backward equation ([31], Section 8.1). In this case, an expected 
value representation in the presence of boundaries is well known. More gen- 
erally, the Feynman-Kac ([31], Section 8.2) formula, at least for linear equa- 
tions with a potential term, has been successfully used in this situation. 
A certain version of this method (Section 3.1), without making the usual 
time reversal substitution, is essentially the same as the method of random 
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At{x,u}i) At{x, UJ2) At{x,u}3) A„,^t{x,uii) 



Fig. 1. Three sample realizations of A without boundaries (left) and with boundaries 
(right). 

characteristics. It is this version that will yield the natural generalization 
of (2.3)-(2.5) in domains with boundary (Theorem 3.1). Before turning to 
the Navier-Stokes equations, we provide a brief discussion on the relation be- 
tween the Feynman-Kac formula and the method of random characteristics. 

3.1. The Feynman-Kac formula and the method of random characteris- 
tics. Both the Feynman-Kac formula and the method of random character- 
istics have their own advantages and disadvantages: The method of random 
characteristics only involves forward SDE's and obtains the solution of (2.8) 
at time t with only the knowledge of the initial data and "X at time t" (or 
more precisely, the solution at time t of the equation (2.3) with initial data 
specified at time 0). However, this method involves computing the spatial 
inverse of X, which analytically and numerically involves an additional step. 

On the other hand, to compute the solution of (2.8) at time t via the 
probabilistic representation using the Kolmogorov backward equation (or 
equivalently, the Feynman-Kac formula with a potential term) when u 
is time dependent involves backward SDE's and further requires the knowl- 
edge of the solution to (2.3) with initial conditions specified at all times 
s <t. However, this does not require computation of spatial inverses and, 
more importantly, yields the correct formulation in the presence of spatial 
boundaries. 

Now, to see the relation between the method of random characteristics 
and the Feynman-Kac formula, we rewrite (2.3) in integral form and keep 
track of solutions starting at all times s > 0. For any s > 0, we define the 
process {Xs^t}t>s to be the flow defined by 

(3.1) Xs,t{x) = X + j^UrO Xs,r{x) dr + V2^{Wt - Ws). 

Now, as always, we let As^t=X~l . Then formally composing (3.1) with Ag^t 
and using the semigroup property X^^t ° X^-^s = X^^t gives the self-contained 
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backward equation for Ag 



(3.2) As,t{x) = x- / UroAr,t{x)dr-V2u{Wt-Ws). 

J s 

Now (2.7) can be written as 

(3.3) et = 9oo Ao,t 

and using the semigroup property A^^s ° ^s,t = ^r,t we see that 

(3.4) et = eso As,t. 

This formal calculation leads to a natural generalization of (2.7) in the 
presence of boundaries. As before, let D C M'^ be a domain with Lipschitz 
boundary and assume, for now, that u is a Lipschitz function defined on all 
of W^. Let Ag^t be the flow defined by (3.2) and for x G D, we define the 
backward exit time (Tt{x) by 

(3.5) at{x) = ml{s\s G [0,t] and Vr G (s, t], A^^f (x) G D]. 

Let g : dD x [0, oo) — )• M and OQ-.D^Mhe two given (regular enough) func- 
tions and define the process 0t by 

9at{x)° ^at{x)A^)^ ifcrt(x)>0, 

%oAo^t{x), iffTt(x) = 0. 

Note that when cTt(x) > 0, equation (3.6) is consistent with (3.4). Thus, (3.6) 
is the natural generalization of (2.7) in the presence of spatial boundaries 
and we expect Of = E9t satisfies the PDE (2.8) with initial data = 
and boundary conditions 9 = g on dD x [0,oo). Indeed, this is essentially 
the expected value representation obtained via the Kolmogorov backward 
equations. 

If an extra term ct{x)9t{x) is desired on the left-hand side of (2.8), then 
we only need to replace (3.6) by 

exp I - / Cs{As,t) dsj g^^(^^) o A^^(^^)^t{x), if <tj(x) > 0, 

V Jcrt(x) J 



(3.6) 9t{x) 



t{x) 



exp 



Jat{x) J 

j Cs{As^t)ds^9oo Ao,t{x), if(Tt(x)=0 



provided c is bounded below. This is essentially the Feynman-Kac formula 
and its application to the Navier-Stokes equations is developed in the next 
section. 

Note that the backward exit time a is usually discontinuous in the spatial 
variable. Thus, even with smooth g,9Q, the process 9 need not be spatially 
continuous. As mentioned earlier, equation (2.6) will now become mean- 
ingless and we will not be able to obtain a SPDE for 0. However, equa- 
tion (2.8), which describes the evolution of the expected value 9 = E9, can 
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be directly derived using the backward Markov property and Ito's formula 
(see, e.g., [18]). We will not provide this proof here but will instead provide 
a proof for the more complicated analogue for the Navier-Stokes equations 
described subsequently. 

3.2. Application to the Navier-Stokes equations in domains with bound- 
ary. First note that if g = in (3.6), then the solution to (2.8) with initial 
data ^0 and 0-Dirichlet boundary conditions will be given by 



(3.7) Ot = Ex{at=o}0o ° Ao^t [i.e., et{x)=Ex{at{x)=o}%°Ao,t{x)]. 



Recall the no-slip boundary condition for the Navier-Stokes equations 
is exactly a 0-Dirichlet boundary condition on the velocity field. Let u be 
a solution to the Navier-Stokes equations in D with initial data uq and 
no-slip boundary conditions. Now, following (3.7), we would expect that 
analogous to (2.5), the velocity field u can be recovered from the flow A^^t 
[equation (3.2)], the backward exit time at [equation (3.5)] and the initial 
data Uq by 



This, however, is false. In fact, there are two elementary reasons one 
should expect (3.8) to be false. First, absorbing Brownian motion at the 
boundaries will certainly violate incompressibility. The second and more 
fundamental reason is that experiments and physical considerations lead us 
to expect production of vorticity at the boundary. This is exactly what is 
missing from (3.8). The correct representation is provided in the following 
result. 

Theorem 3.1. Let u £ C\[0,T);C^{D)) n C{[0,T];C^{D)) be a solu- 
tion of the Navier-Stokes equations (1.1) and (1.2) with initial data uq and 
no-slip boundary conditions. Let A be the solution to the backward SDE (3.2) 
and a be the backward exit time defined by (3.5). There exists a function 
w:dD X [0,T]-^M.^ such that for 



Conversely, given a function w:dD x [0,T] — )• R'^, suppose there exists 
a solution to the stochastic system (3.2), (3.9), (3.10). If further u € C^([0,T); 
C'^{D)) n C([0, T]; C^(Z))), then u satisfies the Navier-Stokes equa- 
tions (1.1)-(1.2) with initial data uq and vorticity boundary conditions 



(3.8) 




when at = 0, 
when fJt > 0, 



we have 



(3.10) 



Ut = PEwt- 



(3.11) 



V X n = V X Ew 



on dD X [0,r]. 
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The proof of Theorem 3.1 is presented in Section 5. We conclude this 
section with a few remarks. 

Remark 3.2. BjV*A^^i^^)^t{x) inequation (3.9) we mean ^*As^t{x)]s=at{x) 
That is, V*^o-t(2:),t(2^) refers to the transpose of the Jacobian of A, evaluated 
at initial time (Tt(x), final time t and position x (see [22, 23, 25] for existence). 
This is different from the transpose of the Jacobian of the function ^crt( ),t(") 
which does not exist as the function is certainly not differentiable in space. 

Remark 3.3 (Regularity assumptions). In order to simplify the pre- 
sentation, our regularity assumptions on u are somewhat generous. Our 
assumptions on u will immediately guarantee that u has a Lipschitz exten- 
sion to M*^. Now the process A, defined to be a solution to (3.2) with this 
Lipschitz extension of u, can be chosen to be a (backward) stochastic flow 
of diffeomorphisms [25]. Thus, is well defined and further defining a 
by (3.5) is valid. Finally, since the statement of Theorem 3.1 only uses val- 
ues of Ag^t for s > at, the choice of the Lipschitz extension of u will not 
matter. See also Remark 5.3. 

Remark 3.4. Note that our statement of the converse above does not 
explicitly give any information on the Dirichlet boundary values of u. Of 
course, the normal component of u must vanish at the boundary of D since u 
is the Leray-Hodge projection of a function. But an explicit local relation 
between w and the boundary values of the tangential component of u can- 
not be established. We remark, however, that while the vorticity boundary 
condition (3.11) is somewhat artificial, it is enough to guarantee uniqueness 
of solutions to the initial value problem for the Navier-Stokes equations. 

Remark 3.5 (Choice of w). We explain how w can be chosen to ob- 
tain the no-slip boundary conditions. We will show (Lemma 5.1) that for w 

defined by (3.9), the expected value w *== Ew solves the PDE 

(3.12) dtwt + {ut ■ V)wt - v^wt + {V*ut)wt = 
with initial data 

(3.13) w\t=o = UQ. 

As shown before, V*ut in (3.12) denotes the transpose of the Jacobian of ut- 
Now, \iu = P?D, then we will have Vxu = Vxu;inD and by continuity, on 
the boundary of D. Thus, to prove existence of the function w, we solve the 
PDE (3.12) with initial conditions (3.13) and vorticity boundary conditions 

(3.14) V xwt = V -K ut on dD. 

We chose w to be the Dirichlet boundary values of this solution. 
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To elaborate on Remark 3.5, we trace through the influence of the vorticity 
on the boundary on the velocity in the interior. First, the vorticity at the 
boundary influences w by entering as a boundary condition on the first 
derivatives for the PDE (3.12). Now, to obtain u we need to find w, the 
(Dirichlet) boundary values of (3.12) and use this to weight trajectories that 
start on the boundary of D. The process of finding w is essentially passing 
from Neumann boundary values of a PDE to the Dirichlet boundary values 
which is usually a nonlocal pseudo-differential operator. Thus, while the 
procedure above is explicit enough, the boundary vorticity influences the 
interior velocity in a highly implicit, nonlocal manner. 

Remark 3.6 (Uniqueness of w). Our choice of w is not unique. Indeed, 
if and w)^ are two solutions of (3.12)-(3.14), then we must have — lu'^ = 
Vq, where q satisfies the equation 

(3.15) V{dtq+{u-V)q-i^Aq) = 

with initial data Vgo = 0. Since we do not have boundary conditions on q, 
we can certainly have nontrivial solutions to this equation. Thus, our choice 
of w is only unique up to addition by the gradient of a solution to (3.15). 

4. Backward Ito integrals. While the formulation of Theorem 3.1 in- 
volves only regular (forward) Ito integrals, the proof requires backward Ito 
integrals and processes adapted to a two parameter filtration. The need 
for backward Ito integrals stems from equation (3.2) which, as mentioned 
earlier, is the evolution of A, backward in time. This is, however, obscured 
because our diffusion coefficient is constant making the martingale term ex- 
actly the increment of the Wiener process and can be explicitly computed 
without any backward (or even forward) Ito integrals. 

To elucidate matters, consider the fiow X' given by 

(4.1) X'^t{a) = a+ j UroX',,.{a)dr+ j ar o X'^^,.{a) dWr- 

J s J s 

If, as usual, A'^^ = (X^^)"^, then substituting formally^ a = A'gf{x) and 
assuming the semigroup property gives the equation 

(4.2) A'^t{x)=x- f UroA'^t{x)dr- f ar o A'^.t{x) dWr 

J s J s 

for the process A'^ The need for backward Ito integrals is now evident; the 
last term above does not make sense as a forward Ito integral since A[.^ is 



*The formal substitution does not give the correct answer when a is not spatially 
constant. This is explained subsequently and the correct equation is (4.3) below. 
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not J> measurable. This term, however, is well defined as a backward Ito 
integral; an integral with respect to a decreasing filtration where processes 
are sampled at the right endpoint. Since forward Ito integrals are more pre- 
dominant in the literature, we recollect a few standard facts about backward 
Ito integrals in this section. A more detailed account, with proofs, can be 
found in [18, 25], for instance. 

Let {Q,T,P) be a probability space, {VFt}f>o be a d-dimensional Wiener 
process on Q and let Ts,t be the cr-algebra generated by the increments Wf — 
Wg' for all s < s' <t' <t, augmented so that the filtration {J^s,t}o<s<t satis- 
fies the usual conditions? Note that for s < s' <t' <t, we have J^s',t' C Ts,t- 
Also Wt — Ws is J%j^t-measurable and is independent of both the past J^o,s, 
and the future Tt,oo- 

We define a (two parameter) family of random variables {Cs,i}o<s<i to 
be a (two parameter) process adapted to the (two parameter) filtration 
{J^s,t}o<s<t, if for all < s < i, the random variable S,s,t is -7^5,rmeasurable. 
For example, S,s,t = Wt — Wg is an adapted process. More generally, if u and a 
are regular enough deterministic functions, then the solution {X'^ t}o<s<t of 
the (forward) SDE (4.1) is an adapted process. 

Given an adapted (two parameter) process and any t > 0, we define the 
backward Ito integral J (,r,tdWr by 



where P = {r = to < ti ■ ■ ■ < tN = t) is a partition of [r,t] and ||P|| is the 
length of the largest subinterval of P. The limit is taken in the sense, 
exactly as with forward Ito integrals (see, e.g., [21], page 148, [27], page 35, 
[25], page 111). 

The standard properties (existence, Ito isometry, martingale properties) 
of the backward Ito integral are, of course, identical to those of the forward 
integral. The only difference is in the sign of the Ito correction. Explic- 
itly, consider the process {A'^ j }o<s<t satisfying the backward Ito differential 
equation (4.2). If {fs,t}o<s<t is adapted, in space and continuously dif- 
ferentiable with respect to s, then the process Bg^t = fs,t ° satisfies the 
backward Ito differential equation 



^By "usual conditions" in this context, we mean that for all s > 0, J-s,s contains all 
^0,00-null sets. Further, Ts,t is right-continuous in t and left-continuous in s. See [21], 
Definition 2.25, for instance. 





* r 
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where = a'^^al^ with the Einstein sum convention. 

Though we only consider solutions to (4.1) for constant diffusion coeffi- 
cient, we briefly address one issue when a is not constant. Our motivation 
for the equation (4.2) was to make the substitution x = A'^ j(x) and formally 
use the semigroup property. This, however, does not yield the correct equa- 
tion when a is not constant and the equation for A'^i = (X^^)"^ involves 
an additional correction term. To see this, we discretize the forward integral 
in (4.1) (in time) and substitute a = A'^ i{x). This yields a sum sampled at 
the left endpoint of each time step. While this causes no difficulty for the 
bounded variation terms, the martingale term is a discrete approximation 
to a backward integral and hence, must be sampled at the right endpoint of 
each time step. Converting this to sum sampled at the right endpoint via 
a Taylor expansion of a is what gives this extra correction. Carrying through 
this computation (see, e.g., [25], Section 4.2) yields the equation 

A'^i-{x) = x- f Ur o A'j. ^{x) dr - f ar o A'^ ^{x) dWr 

(4.3) 

+ / {djai''' o A'^ t{x)){4''' o A'^ t{x))eidr, 



where {ej}i<j<d are the elementary basis vectors and a^'^ denotes the i,jth 
entry in the d x d matrix a. 

We recall that the proof of the (forward) Ito formula involves approxi- 
mating / by its Taylor polynomial about the left endpoint of the partition 
intervals. Analogously, the backward Ito formula involves approximating / 
by Taylor polynomial about the right endpoint of partition intervals, which 
accounts for the reversed sign in the Ito correction. 

Finally, we remark that for any fixed t > 0, the solution {As^t}o<s<t of the 
backward SDE (3.2) is a backward strong Markov process [the same is true 
for solutions to (4.3)]. The backward Markov property states that r < s <t 
then 

Ers,tf° Ar,t{x) = EA,^t(x)f O AA^) = [Ef ° A,s(?/)]y=A.,t{x), 

where Ejr^ ^ denotes the conditional expectation with respect to the c-alge- 
bra Ts,t and E^^ t(x) the conditional expectation with respect to the cr-alge- 
bra generated by the process Ag^tix)- 

For the strong Markov property (we define cr to be a backward t -stopping 
time^ if almost surely a < t) and for all s <t, the event {a > s} is J-g^t- 
measurable. Now if a is any backward t-stopping time with r < a <t almost 



®Our use of the term backward f-stopping time is analogous to s-stopping time in [18], 
page 24. 
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surely, the backward strong Markov property states 

Et.J° ArA^) = Ea^J o Ar,t{x) = [Ef o Ar,s{y)] s=a, 

y=A^^t{x). 

The proofs of the backward Markov properties is analogous to the proof of 
the forward Markov properties and we refer the reader to [18], for instance. 

5. The no-slip boundary condition. In this section we prove Theorem 3.1. 
First, we know from [22, 23] that spatial derivatives of A can be interpreted 
as the limit (in probability) of the usual difference quotient. In fact, for 
regular enough velocity fields u (extended to all of M"'), the process A can, 
in fact, be chosen to be a flow of diffeomorphisms of M'^ (see, e.g., [25]) 
in which case A is surely differentiable in space. Interpreting the Jacobian 
of A as either the limit (in probability) of the usual difference quotient or as 
the Jacobian of the stochastic flow of diffeomorphism, we know [22, 23, 25] 
that VA satisfies the equation 

(5.1) VAsAx) = I- [ Vur\ArA^)yA,t{x)dr, 

J s 

obtained by formally differentiating (3.2) in space. Here I denotes the d x d 
identity matrix. We reiterate that equation (5.1) is an ODE as the Wiener 
process is independent of the spatial parameter. 

Lemma 5.1. Let D,u,T be as in Theorem 3.1, a he the backward exit 
time from D [equation (3.5)] and A he the solution to (3.2) with respect to 
the backward stopping time a. 

(1) Let -u) G Ci([0,r);C2(D)) n C([0,r];Ci(i5)) be the solution of (3.12) 
with initial data (3.13) and boundary conditions 

(5.2) w = w on dD. 

Then, for w defined by (3.9), we have w = Ew. 

(2) Let w he defined by (3.9) and w = Ew as above. If for all t £ (0,T], 
idt € T>(A.A CLnd w is in time, then w satisfies 

(5.3) dtw + Ltw + {V*u)w = Q, 
where Lt is defined by 

(5.4) = hm ^(^)-^^(-4w..(.),(x)) 

s^t- t — S 

and T>{A.A is the set of all (j) for which the limit on the right-hand side 
exists. Further, w has initial data uq and boundary conditions (5.2). 
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Before proceeding any further, we first address the relationship between 
the two assertions of the lemma. We claim that if w ^ C^((0, T); C^(L')), 
then equation (5.3) reduces to equation (3.12). This follows immediately 
from the next proposition. 

Proposition 5.2. //(/> e C^(Z)), then for any t € {0,T], (f) £'D{A.^t) and 
further, Lfcj) = {uf V)(j) — i^Acp. 

Proof. Omitting the spatial variable for notational convenience, the 
backward ltd formula gives 

(f>- (f>0 As\/at,t = (P° At^t -(t>° As\/at,t 



Since s V fj^ is a backward t-stopping time, the second term above is a mar- 
tingale. Thus 



since the process A has continuous paths and at <t on the interior of D. 



Remark 5.3. One can weaken the regularity assumptions on u in the 
statement of Theorem 3.1 by instead assuming for all t G (0,T], ut £ "D^A.^t) 
and is in time, as with the second assertion of Lemma 5.1. However, while 
the formal calculus remains essentially unchanged, there are a couple of tech- 
nical points that require attention. First, when assumptions on smoothness 
of u up to the boundary is relaxed (or when dD is irregular), a Lipschitz 
extension of u need not exist. In this case, we can no longer use (3.5) to 
define a. Further, we can not regard the process yl as a stochastic flow 
of diffeomorphisms and some care has to be taken when differentiating it. 
These issues can be addressed using relatively standard techniques and once 
they are sorted out, the proof of Theorem 3.1 remains unchanged. 

Now we prove the first assertion of Lemma 5.1. 

Proof. Recall that \/*As^t is differentiable in s. Differentiating (5.1) in s 
and transposing the matrices gives 





{ut • V)0 - uA(j) 



□ 



(5.5) 



dsVAs,t{x) = V*As,t{x)V*Us\AsA^). 
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Let t £ {0,T], X £ D and a' be any backward t-stopping time with a' > at{x) 
almost surely. Omitting the spatial variable for convenience, the backward 
Ito formula and equations (3.12) and (5.5) give 

= V*At^t'Wt O At^t - V*A^,^tWa' O Aa',t 

= / drV* Ar^tWr O Ar^t 
Ja' 

+ / V* Ar^t{drWr + {Ur ■ V)Wr — f/^Wr) o Ar^t dv 
Ja' 

+ V2^ [ {V*Art){V*Wr)oArtdWr 
Ja' 

= / V* Ar^t{{^*Ur)Wr + drWr + {Ur ■ V)Wr — l^i^Wr) o Ar,tdr 
Ja' 

+ V2^ [ {V*Ar,t){V*Wr)oAr,tdWr 
Ja' 

= V2u [ {V*Ar,t)Wr o Ar,t dWr- 
Ja' 

Thus, taking expected values gives 

(5.6) wt{x) = EV*A^,^t{x)wa' o A^,^t{x). 

Recall that when at{x) > 0, Ao-j(^) j(x) G dD. Thus, choosing a' = at{x) and 
using the boundary conditions (5.2) and initial data (3.13), we have 

(5.7) - I ^ ^^^^^^^^^ ^ Q_ 
Substituting this in (5.6) completes the proof. □ 

In order to prove the second assertion in Lemma 5.1, we will directly 
prove (5.6) using the backward strong Markov property. Before beginning 
the proof, we establish a few preliminaries. 

Let u, T, a, A, w, w be as in the second assertion of Lemma 5.1. Given 
X £ D and a dx d matrix M, define the process {Bg^tix, M)}^^(^^-^.cs<t<T to 
be the solution of the ODE 

Bs,t{x,M) = M- I Vur\ArA^)BrA^,M)dr. 

J s 

If / denotes the dx d identity matrix, then by (5.1) we have Bs^t{x,I) = 
V As^t{.x) for any (Jt{x) < s <t <T. Further, since the evolution equation 
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Bs,t{x, M) = Bs,t{xJ)M = \IAs,t{x)M. 



Note that for any fixed t G (0,T], tlie process As^t}o<s<t is not a backward 
Markov process. Indeed, the evolution of V Ag^t at any time s <t depends on 
the time s through the process Ag^t appearing on the right-hand side in (5.1). 
However, process {As^t-,'^ ■^s,t) [or equivalently the process (^s^t, is 
a backward Markov process since the evolution of this system now only 
depends on the state. This leads us to the following identity which is the 
essence of proof of the second assertion in Lemma 5.1. 

Lemma 5.4. Choose any backward t-stopping time a' with a' > at{x) 
almost surely. Then 



holds almost surely. 

This follows from an appropriate application of the backward strong 
Markov property. While this is easily believed, checking that the strong 
Markov property applies in this situation requires a little work and will 
distract from the heart of the matter. Thus, we momentarily postpone the 
proof of Lemma 5.4 and proceed with the proof of the second assertion of 
Lemma 5.1. 

Proof of Lemma 5.1. We recall w = Ew where w is defined by (3.9). 
By our assumption on u and dD, the boundary conditions (5.2) and initial 
data (3.13) are satisfied. For convenience, when y G dD, t > 0, we define 
wt{y) = w{y) and when t = 0, y £ D, we define WQ{y) = ito(y)- 

Let X £ D, t £ (0,T] as used before. Let a' be any backward t-stopping 
time with a' > (Tt{x). First, if a' = at{x) almost surely, then, since the 

point {A^^(^ri.^^i,t) belongs to the parabolic boundary dp{D x [0,T]) {dD x 
[0,r]) U {D X {0}), our boundary conditions and initial data will guaran- 



Now, for arbitrary a' > (Tt{x), we will use Lemma 5.4 to deduce (5.6) 
directly. Indeed, 



(5.9) 



Et^,^, Bl^ (x, I)w„^ (^) o A„^ (^)^j (x) 



M=B^,,{x,I) 



tee (5.6). 



wt{x) = £;VM^^(^)_t(x)'u)^^(^) o 

= EE^^,^^ Kt{x),t i^^ I)Wat (x) ° Kt {x),t (X) 
= ^([^^:.fe),r(y>^^)^<x.fe)0^a.fe),r(?/)] 



\r=cT',y=A^,^[x), 
M=B^,,{x,I) 
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M=B^,^,{x,I) 

showing that (5.6) holds for any backward t stopping time a' > at{x). 

Now, choose a' = s V at{x) for s <t. Note that for any x € D, we must 
have at{x) < t almost surely. Thus, omitting the spatial coordinate for con- 
venience, we have 

= lim ^ — — = lim —^{wt - EV* As\jat,tWs\jat ° As\jat,t) 

s^t- t — S s^t- t — S 

= lim (—^[wt - Ewt o Ass/„ t] 
s^t- \t — s 

+ —^E{wt - WsS/at) ° Assjat,t 

t — s 

= Ltwt + dtwt + {\/*ut)wt, 
on the interior of D. The proof is complete. □ 

It remains to prove Lemma 5.4. 

Proof of Lemma 5.4. Define the stopped processes ^s^((x)=j4o-((^.)vs,t(2;) 
and B'g ^{x,M) = B^^(^^-^-^g^-i-{x, M). Define the process C by 

Cs,t{x, M, t) = {A'^/x),B'^^t{x, M),T + t- at{x) V s). 

Note that for any given s < t, we know that at{x) need not be J^g^t measur- 
able. However, at{x) V s is an J^g^t measurable backward t-stopping time. 
Thus, A'g^, B'gi and, consequently, Cg^t are all Fg^t measurable. 

Now we claim that almost surely, for 0<r<s<t<T, we have the 
backward semigroup identity 

(5.10) Cr,gOCg^t = Cr,t. 

To prove this, consider first the third component of the left-hand side of (5.10): 

(5.11) eg) o Cg,t{x, M,T) = {T + t- atix) V s) + s - a,(A', ^(x)) V s. 

Consider the event {s > at{x)}. By the semigroup property for A and strong 
existence and uniqueness of solutions to (3.2), we have ag{Ag^t{x)) = at{x) 
almost surely. Thus, almost surely on {s > at{x)}, we have 

eg) o Cg^t{x, M,T) = {T + t-s) + s- at{x) V s 

= T + t-at{x)Vr = eg) (x, M, r) . 
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Now consider the event {s<at}. Weknow A'gj-{x)£dD and so cTs(A'^ j(x)) = s. 
This gives 

eg o Cs,tix, M,T) = {T + t- atix)) + s-s = T + t- at{x) V r = cij(x) 

almost surely on {s < at{x)}. Therefore, we have proved almost sure equality 
of the third components in equation (5.10). 

For the first component cf^"* = A'^ j, consider as before the case s > (Tt{x). 
In this case A'^ ^ = As^t and the semigroup property of A gives equality of the 
first components in (5.10) almost surely on {s > at{x)}. When s < (7t{x), as 
before, A'^ ^^ € dD and c7s(A'^ j(x)) = s. Thus, 

almost surely on s < at{x). This shows almost sure equality of the first com- 
ponents in equation (5.10). Almost sure equality of the second components 
follows similarly, completing the proof of (5.10). 

Now, for < r < s < t <T, the random variable Cg^t is J^s,t measurable 
and so must be independent of Jv^g. This, along with (5.10), will immediately 
guarantee the Markov property for C. Since the filtration J^.^. satisfies the 
usual conditions and for any fixed t the function s i— )• Cs,t is continuous, C 
satisfies the strong Markov property (see, e.g., [18], Theorem 2.4). 

Thus, for any fixed t G [0, T] and any Borel function ip, the strong Markov 
property gives 

E^^,^ip{Co,t{xJ,0)) = [Eip{Cr,tiy,M,T))]r=a', 

{y,M,r)=Co^,,{x,I,0) 

= [E(p{Cr,t{y,M,T))] 

M=B^, ((x,/),r=(7r{x), 

almost surely for any xGM'^, MsR'^^, r>0. Choosing ip{x, M, T)=M*wt-T{x) 
proves (5.9). □ 

Now a direct computation shows that if w satisfies (3.12), then u = Pw) 
satisfies (1.1) regardless of our choice of w. Of course, we will only get the 
no-slip boundary conditions with the correct choice of w. We first obtain 
the PDE for u. 

Lemma 5.5. If w satisfies (3.12) and u = Fw, then u satisfies (1.1) 
and (1.2). 

Proof. By definition of the Leray-Hodge projection, u = w + Vq for 
some function q and equation (1.2) is automatically satisfied. Thus, using 



STOCHASTIC-LAGRANGIAN NAVIER-STOKES WITH BOUNDARY 



19 




Vp = V(i|n|2 + dtqt + {ut ■ V)qt - uAqt), 



equation (5.12) becomes (1.1). □ 

Now to address the no-slip boundary condition. The curl of w satisfies the 
vorticity equation which is how the vorticity enters our boundary condition. 

Lemma 5.6. Let w be a solution of (3.12). Then ^ = V x tj) satisfies the 
vorticity equation 



Proof. We only provide the proof for d = 2>. For this proof we will use 
subscripts to indicate the component instead of time as we usually do. If 
/c € {1,2,3} are all distinct, let £ijk denote the signature of the permu- 
tation (1, 2, 3) I— 7- (i, j, k). For convenience, we let Eijk = if i, j, k are not all 
distinct. Using the Einstein summation convention, ^ = V x "iD translates to 

= ^ijkdjWk on components. Thus, taking the curl of (3.12) gives 

(5.14) dtCi + (u ■ V)^i - l/A^i + Eijk djUm dmWk + £ijk dkUm djWm = 

because £ijkdjdkUmWm = 0. Making the substitutions j ^k and k i-^ j in 
the last sum above we have 

dmWk + £ijk dkUm djWm = Sijk djUm{dmWk - dkWm) 



where Sij denotes the Kronecker delta function and the last equality follows 
because djuj =0. Thus, (5.14) reduces to (5.13). □ 

Theorem 3.1 now follows from the above lemmas. 

Proof of Theorem 3.1. First, suppose u is a solution of the Navier- 
Stokes equations, as in the statement of the theorem. We choose w as ex- 
plained in Remark 3.5. Notice that our assumptions on u and D will guar- 
antee a classical solution to (3.12)-(3.14) exists on the interval [0,r] and 
thus, such a choice is possible. 



(5.13) 




ifd = 2, 
ifd = 3. 




,n 
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By Lemma 5.1 we see that for w defined by (3.9), the expected value w = 
Ew satisfies (3.12) with initial data (3.13) and boundary conditions (5.2). 
By our choice of w and uniqueness to the Dirichlet problem (3.12), (3.13) 
and (5.2), we must have the vorticity boundary condition (3.14). 

Now, let ^ = V X "iD and w = V x u. By Lemma 5.6, we see that ^ satisfies 
the vorticity equation (5.13). Since u satisfies (1.1) and (1.2), it is well known 
(see, e.g., [15, 28] or the proof of Lemma 5.6) that oj also satisfies 



From (3.14) we know ^ = a; on dD x [0,T]. By (3.13), we see that .^o = 
V X = Wo and hence, ^ = uj on the parabolic boundary dp{D x [0, T]). 

The above shows that uj and both satisfy the same PDE [equations (5.13) 
or (5.15)] with the same initial data and boundary conditions and so we must 
have ^ = u; on D X [0, T] . Thus, V x ?D = V x n in D x [0, T] showing u and w 
differ by a gradient. Since V • n = and u = on dD x [0, T] , we must have 
u = PtD proving (3.10). 

Conversely, assume we have a solution to the system (3.2), (3.9) and (3.10). 
As stated above, Lemma 5.1 shows w = Ew satisfies (3.12) with initial 
data (3.13). By Lemma 5.5 we know u satisfies the equation (1.1) and (1.2) 
with initial data uq. Finally, since equation (3.10) shows Vxii = Vxu;in 
D X [0,T] and by continuity, we have the boundary condition (3.11). □ 

6. Vorticity transport and ideally conserved quantities. The vorticity is 
a quantity which is of fundamental importance, both for the physical and 
theoretical aspects of fluid dynamics. To single out one among the numerous 
applications of vorticity, we refer the reader to two classical criterion which 
guarantee global and existence and regularity of the Navier-Stokes equations 
provided the vorticity is appropriately controlled: the first due to Beale, Kato 
and Majda [2] and the second due to Constantin and Fefferman [9]. 

For the Euler equations, exact identities and conservation laws governing 
the evolution of vorticity are well known. For instance, vorticity transport 
[equation (6.1)] shows that the vorticity at time t followed along streamlines 
is exactly the initial vorticity stretched by the Jacobian of the flow map. 
Similarly, the conservation of circulation [equation (6.9)] shows that the line 
integral of the velocity (which, by Stokes theorem, is a surface integral of 
the vorticity) computed along a closed curve that is transported by the fluid 
flow is constant in time. 

Prior to [11], these identities were unavailable for the Navier-Stokes equa- 
tions. In [11], the authors provide analogues of these identities for the Navier- 
Stokes equations in the absence of boundaries. These identities, however, do 
not always prevail in the presence of boundaries. 



(5.15) 




if d = 2, 
if d = 3. 



STOCHASTIC-LAGRANGIAN NAVIER-STOKES WITH BOUNDARY 



21 



In this section we illustrate the issues involved by considering three in- 
viscid identities. All three identities generalize perfectly to the viscous sit- 
uations without boundaries. In the presence of boundaries, the first iden- 
tity (vorticity transport) generalizes perfectly, the second identity (Ertel's 
Theorem) generalizes somewhat unsatisfactorily and the third identity (con- 
servation of circulation) has no nontrivial generalization in the presence of 
boundaries. 

6.1. Vorticity transport. Let be a solution to the Euler equations with 
initial data uq. Let the inviscid flow map defined by (2.1) and for any 
t > 0, let = {X^)~^ be the spatial inverse of the diffeomorphism . The 
vorticity transport (or Cauchy formula) states 



where we recall that the vorticity a; is defined hy lo —V x u and where 
ujq = V X Uq is the initial vorticity. 

In [11], the authors obtained a natural generalization of (6.1) for the 
Navier-Stokes equations in the absence of spatial boundaries. If u solves (1.1) 
and (1.2) with initial data uq and X is the noisy flow map defied by (2.3)- 
(2.4), then a; = V x u is given by 



We now provide the generalization of this in the presence of boundaries. 
Note that for any t>0, (VXt) oAt = (VAt)'^, so we can rewrite (6.2) com- 
pletely in terms of the process A. Now, as usual, we replace A = X~^ with 
the solution of (3.2) with respect to the minimal existence time a. We recall 
that in Theorem 3.1, in addition to "starting trajectories at the boundary," 
we had to correct the expression for the velocity by the boundary values 
of a related quantity (the vorticity). For the vorticity, however, we need 
no additional correction and the interior vorticity is completely determined 
given A, a and the vorticity on the parabolic boundary^ dp{D x [0,T]). 

Proposition 6.1. Letu be a solution to (1.1) and (1.2) in D with initial 
data Uq and suppose w = V x n G C^{[^,T);C'^{D)) n C([0,T] x D). Let uj 
denote the values of uj on the parabolic boundary dp{D x [0,T]). Explicitly, 
u is defined by 



(6.1) 




if (i = 2, 
if d = 3, 



(6.2) 




if d = 3. 




if X & D and t = 
if X ^ dD . 



■'Recall the parabolic boundary dp{D x [0, T]) is defined to be {D x {0}) U {dD x [0, T)). 
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Then. 



(6.3) ut{x) = 



(^<7t(x),t(2;))] 



tfd = 2, 
ifd = 3. 



Proposition 6.2. More generally, suppose tD is any function defined on 
the parabolic boundary of D x [0,T] and let uj be defined by (6.3). If for all 
t G (0,T], ujt G ^{A.^t) CLnd UJ is in time, then to satisfies 



with u = Cj on the parabolic boundary. Here, Lt is the generator of A.^t! 
'D(A.^t) is the domain of Lt. These are defined in the statement of Lemma 5.1. 

Of course, Proposition 6.2, along with Proposition 5.2 and uniqueness 
of (strong) solutions to (5.15), will prove Proposition 6.1. However, direct 
proofs of both Proposition 6.2 and Proposition 6.1 are short and instructive 
and we provide independent proofs of each. 

Proof of Proposition 6.1. We only provide the proof when ci = 3. 
As shown before, differentiating (5.1) in space and taking the matrix inverse 
of both sides gives 



almost surely. Now choose any x £ D, t > and any backward t-stopping 
time a' > (Tt{x). Omitting the spatial parameter for notational convenience, 
the backward Ito formula gives 




ifd = 2 
ifd = 3 



(6.4) 



dr{VAr,t{x)) 



-1 



LUt - (VA^f^t) uj^, o A^,^t 

= {VAt,ty^ujt o At,t - {VA^>^t)~^iOa' o A^,^t 




+ / (VAr^t) {drUJr + [Ur ' V)^,. — uAuJr) O A^-^t dr 
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Ja' 

= V2u [ {VArty^{VuJr)oArtdWr. 
Ja' 

Thus, taking expected values gives 

(6.5) ujt = E[{VA„,^t)~^uja'°A,,^t]. 

Choosing a' = (Jt{x) and using the fact that Ao-t(x),t(^) always belongs to 
the parabolic boundary finishes the proof. □ 

Proof of Proposition 6.2. Again, we only consider the case d = 3. 
We will prove equation (6.5) directly and then deduce (5.15). Let the pro- 
cess B be as in the proof of the second assertion of Lemma 5.1 and use 
to denote the process consisting of matrix inverses of the process B. Pick 
X & D, t& {0,T] and a backward t-stopping time a' > at{x). Using (5.9) we 
have 

ujtix) = ^[(V^^^(^)^t(x))"^(I'^^(^)(A^^(^,)^t(x))] 

= ee^^,jb;^\,^^/x,i)u^^(^^) o ^<,,(^),t(x)] 

M=B^,,{x,I) 

= E{[M-^EB;^^y^/yJ)u^^^y) o A^^^y)^r{y)]r=a',y=A^,4x),) 

= E[{VA^,^t{x))-^uj^, o A^,^tix)], 
proving (6.5). 

As stated before, choose s <t and a' = at{x) Vs. Omitting the spatial 
parameter for notational convenience gives 

0= lim — — —= lim -^[ujt - EiV A^^ys^tT^^atys ° A„^ys,t] 

s^t- t — S s^t~ t — S 

= lim [ -^[ujt- Eioto Afj^ys^t] 
s^t- \ t — s 

+ E[uJt — (^(TtVs] ° ^o-tVs,t 

t — S 

= LtUJt + dtujt - iVut)ujf □ 
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We remark that the vorticity transport in Propositions 6.1 or 6.2 can be 
used to provide a stochastic representation of the Navier-Stokes equations. 
To see this, first note that the proofs of Propositions 6.1 and 6.2 are inde- 
pendent of Theorem 3.1. Next, since u is divergence free, taking the curl 
twice gives the negative laplacian. Thus, provided boundary conditions on u 
are specified, we can obtain u from uj by 

(6.6) ut = i-Ay^V xut. 

Therefore, in Theorem 3.1 we can replace (3.10) by (6.3) and (6.6), where uj is 
the vorticity on the parabolic boundary and we impose 0-Dirichlet boundary 
conditions on (6.6). 

6.2. Ertel's theorem. As shown above, we use a superscript of to denote 
the appropriate quantities related to the Euler equations. For this section 
we also assume d = 3. Ertel's theorem says that if 9^ is constant along 
trajectories of X^, then so is {uj^ ■ S/)6^. Hence, <j)^ = [uj^ ■ V)6^ satisfies the 
PDE 

at/ + (u ■ V)/ = 0. 

For the Navier-Stokes equations, we first consider the situation without 
boundaries. Let u solve (1.1) and (1.2), X be defined by (2.3), A be the 
spatial inverse of X and define ^ by 

^t{x) = {VAt{x))~^ujooAtix), 

where coq = 'V x uq is the initial vorticity. From (6.2) we know that uj = 
V X li = ES^. Now we can generalize Ertel's theorem as follows: 

Proposition 6.3. Let 6 be a C^(R^) valued process. If 6 is constant 
along trajectories of the (stochastic) flow X , then so is (^ • V)^. Hence, 
4) = E{i- V)9 satisfies the PDE 

(6.7) dt4>t + {ut ■ V)(tJt - i^A(ljt = 0, 
with initial data (cjq • V)^o- 

Proof. If 6 is constant along trajectories of X, we must have 6t = 9o° At 
almost surely. Thus, 

(6 • ^Wt = {V9t)Ct = V9oU {VAt){VAt)-^ujo o At = (Co • V^o) o At, 

which is certainly constant along trajectories of X. The PDE for (p now 
follows immediately. □ 

Now, in the presence of boundaries, this needs further modification. Let A 
be a solution to (3.2) and a be the backward exit time of A from D. The 
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notion of "constant along trajectories" now corresponds to processes 6 de- 
fined by 

(6.8) 9t{x)=0^^^,)iA^^^^)^t), 

for some function 9 defined on the parabolic boundary of D. 

Irrespective of the regularity of D and 0, the process 6 will not be contin- 
uous in space, let alone differentiable. The problem arises because while A is 
regular enough in the spatial variable, the existence time at is not. Thus, we 
are forced to avoid derivatives on a in the statement of the theorem, leading 
to a somewhat unsatisfactory generalization. 

Proposition 6.4. Let 9 he aC^ function defined on the parabolic bound- 
ary of D X [0, T] and let 9' be any extension of 9, defined in a neighborhood 
of the parabolic boundary of D x [0,T]. If 9 is defined by (6.8), then 

Ux) = E[{it ■ V){9i o As,t){x)]s=.,(^,) 

satisfies the PDE (6.7) with initial data (wq • V)^o o-nd boundary conditions 
(j)t{x) = (wt • V)9'{x) for xedD. 

Remark. A satisfactory generalization in the scenario with boundaries 
would be to make sense of E{^t • ^)9t (despite the spatial discontinuity of 9) 
and reformulate Proposition 6.4 accordingly. 

Note that when D = M*^, then cjf = and hence, </)t = E{^t ■ '^)9t- In this 
case Proposition 6.4 reduces to Proposition 6.3. The proof of Proposition 6.4 
is identical to that of Proposition 6.3 and the same argument obtains 

[iCt ■ VM ° ^,t)(x)],=,,(,) = [(6 • ^)0'siy)] s=M-), , 

which immediately implies (6.7). 

6.3. Circulation. The circulation is the line integral of the velocity field 
along a closed curve. For the Euler equations, the circulation along a closed 
curve that is transported by the flow is constant in time. Explicitly, let u^, 
X^, A^, uq be as in the previous subsection. Let F be a rectifiable closed 
curve, then for any t > 0, 

(6.9) * u^fdl= <p ul- dl. 

JxO(r) Jr 

For the Navier-Stokes equations, without boundaries, a generalization of (6.9) 
was considered in [11]. Let u solve (1.1) and (1.2), X be defined by (2.3) 
and (2.4) and A be the spatial inverse of X. Then 

(6.10) (f, ut ■ dl = E (h UQ-dl. 

Jr JAtiv) 
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A proof of this (in the absence of boundaries) follows immediately from 
Theorem 2.1. Indeed, 



where the first equality follows by definition of line integrals, the second 
because the line integral of gradients along closed curves is and the last 
by Fubini and (2.5). 

Equation (6.10) does not make sense in the presence of boundaries, as the 
curves one integrates over will no longer be rectifiable! 
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